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1 Introduction 

Our understanding of the gamma-ray bursts (GRBs) has advanced rapidly 
since the discovery of long-lived "afterglow" emission from these events. Ra- 
dio afterglow studies have become an integral part of this field, providing 
complementary and sometimes unique diagnostics on GRB explosions, their 
progenitors, and their environments. The reason for this is that the radio part 
of the spectrum is phenomenologically rich. This can be illustrated simply by 
calculating the brightness temperature (Tf, oc Y v /{9 s v) 2 ) for a 1 mJy cen- 
timeter wavelength source at cosmological distances (~ 10 28 cm), expanding 
with V exp < c one week after the burst. Since the derived ~ 10 13 K is 
well in excess of the Tic ~ 10 11 — 10 12 K limit imposed by inverse Compton 
cooling, it follows, independent of any specific afterglow model, that the radio 
emission must originate from a compact, synchrotron-emitting source that is 
expanding superluminally (i.e. Tb ~ T x Tic, T >> 1). Likewise, since the 
brightness temperature cannot exceed the mean kinetic energy of the elec- 
trons, the emission is expected to be self-absorbed at longer wavelengths pQ. 
Finally, strong modulation of the centimeter signal is expected on timescales 
of hours and days because the angular size 9 S of this superluminal source is 
comparable to the Fresnel angle of the turbulent ionized gas in our Galaxy 2 . 
Synchrotron self-absorption, interstellar scintillation, forward shocks, reverse 
shocks, jet-breaks, non-relativistic transitions and obscured star formation are 
among the phenomena routinely observed. 

This short review is divided into two parts. The first section (SJ3J) is a sum- 
mary of the current search strategies and the main observational properties of 
radio afterglows. In the second section (@ we highlight the key scientific con- 
tributions made by radio observations, either alone or as part of panchromatic 
studies. By necessity we will restrict this brief review to long-duration GRBs, 
although radio afterglows have also been detected toward the newly classified 
X-ray flashes, and searches have been carried out toward short bursts |3J- 
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2 Detection Statistics and Observational Properties 

The search for a radio afterglow is initiated either by a satellite localization of 
the burst, or by the detection of the X-ray or optical afterglow. The current 
search strategy has been to use the Very Large Array (VLA) 1 or the Australia 
Telescope Compact Array (ATCA; for declinations, 5 < -40°) 2 at 5 GHz or 
8.5 GHz. These frequencies were chosen as a compromise between the need to 
image the typical error box size of 30-100 arcmin 2 , while having the requisite 
sensitivity to detect afterglows at sub-milliJansky levels. At lower frequencies 
the afterglow is attenuated by synchrotron self-absorption (f v oc v 2 ), while at 
higher frequencies the field-of-view is proportionally smaller (FOVoc v 2 ). For 
typical integration times (10 min at the VLA, and 240 min at the ATCA) the 
rms (receiver) noise is 30-50 /xJy. Follow-up observations of detected after- 
glows were carried out by a network of radio facilities at centimeter, millimeter 
and submillimeter wavelengths 0]. 
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Fig. 1. (Left) Histogram distribution of flux densities (or upper limits) at 8.5 GHz 
for a complete sample of bursts. The hatched histogram shows the distribution for 
the detections only. (Right) Histograms of radio luminosity from the same sample 
but restricted to the subset of bursts with known redshifts. The hatched histogram 
shows the distribution for bursts with detected radio afterglows only. 



In the five year period beginning in 1997 and ending in 2001 approximately 
1500 radio flux density measurements (or upper limits) were made toward 75 
bursts |51 ■ From these 75 GRBs, there are a total of 32/36 successful X- 
ray searches, 27/70 successful optical searches, and 25/75 successful searches. 
These afterglow search statistics illustrate a well-known result, namely that 
the detection probability for X-ray afterglows is near unity, while for optical 

The NRAO is a facility of the National Science Foundation operated under co- 
operative agreement by Associated Universities, Inc. 
2 The Australia Telescope is funded by the Commonwealth of Australia for opera- 
tion as a National Facility managed by CSIRO. 
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afterglows and radio afterglows it is 40% and 33%, respectively. The origin 
of these optically "dark bursts" could either be due to intrinsic effects (i.e. 
inadequate search due to rapid evolution of the afterglow and/or an under- 
energetic GRB) [HI El, or an extrinsic effect (i.e. extinction of the optical flux 
caused by circumburst dust or by the intergalactic medium) 

To accurately derive the fraction of "radio quiet" bursts it is necessary to 
incorporate both detections and upper limits in a statistically sound manner. 
This has been done in Fig. where flux density distribution at 8.5 GHz is 
shown for a sample of 44 GRBs, toward which measurements or upper limits 
have been made between 5 and 10 days after a burst. The time since the burst 
is an important variable since radio light curves do not exhibit the simple 
power-law decays seen in X-ray and optical afterglows, but rise to a peak on 
average about one week after the burst and decay on timescales of a month. 
The mean of the 19 detections in Fig. ^is 315±82 /iJy. Adding in the non- 
detections, and using the Kaplan-Meier estimator JOj shifts this to 186±40 
jtxjy. Approximately 50% of all bursts have radio afterglows at 8.5 GHz above 
110 /iJy, while fewer than 10% exceed 500 ^Jy. The relatively small range 
of peak flux densities in Fig. ^ suggests that the fraction of "radio quiet" 
bursts is largely determined by instrumental sensitivities. With the arcsecond 
localizations provided by the Swift satellite (launch in 2004) it will be possible 
to routinely detect all afterglows with centimeter radio emission above 100 
/iJy. Increasing the fraction of detected radio afterglows significantly above 
50% will require the sensitivity improvements provided by the Expanded Very 
Large Array* (complete in 2010). 

From this sample of peak flux densities we also derive the peak spectral 
radio luminosity in Fig. given by L v — AttF^ d\ (1 + z) 1+/3 ~ Q , where F v oc 
t a v" and a = 1/2 and /3 = 1/3 has been assumed, corresponding to an 
optically thin, rising light curve. The GRB redshifts lie in the range between 
z =0.36 to z=4.5. The peak of the distribution is centered on 10 31 erg s^ 1 
Hz -1 and is similar to low-luminosity FRI radio galaxies like M87. More 
interestingly, a comparison between this GRB sample and a sample of Type 
Ib/c supernovae JI] shows that the later is four orders of magnitude less 
luminous. Since radio emission is sensitive to the relativistic energy content of 
the shock, independent of the initial geometry of the explosion, this has been 
used to argue that the majority (<97%) of nearby Type Ib/c supernovae do 
not produce a GRB-like event, such as that seen toward SN 1998bw |12) . 

3 Phenomenology and Interpretation 

In this section we will follow the evolution a GRB and its radio afterglow 
depicted schematically in Fig. The observations span four orders of mag- 
nitude in time (0.1-1000 days) and three orders of magnitude in frequency 

3 http://www.aoc.nrao.edu/evla/ 
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(0.8-660 GHz), so it should be no surprise that radio light curves exhibit a 
rich phenomenology. To interpret these observations we will rely on the highly 
successful "standard fireball model" UIJ. In this model there is an impulsive 
release of kinetic energy (~10 51 erg) from the GRB event which drives an 
ultra-relativistic outflow into the surrounding medium whose hydrodynami- 
cal evolution is governed by the kinetic energy released, the density structure 
of the circumburst medium and the geometry of the outflow. Synchrotron 
emission is produced by this relativistic shock which accelerates electrons to 
a power-law distribution. It is through the study of temporal (and spectral) 
evolution of afterglow light curves that we can gain insight into the physical 
conditions of the shock and the central engine that produced it. 
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Fig. 2. A schematic radio afterglow light curve. Timescales and scalings for the 
temporal evolution are indicated. The list summarizes aspects of the flux evolution 
which are unique to the radio bands (Lorentz factor, F; source size, 9; energy, E; 
density, n; jet opening angle, #j e t; density profile; magnetic field strength, B; and 
obscured star formation rate). 



Despite response times as short as 2 hrs, centimeter searches (§2} are rarely 
successful until a day or more after a burst. Broadband afterglow spectra 
show that centimeter emission is attenuated as a result of synchrotron self- 
absorption ^1]. Typical observed values for the self- absorption frequency v a 
are 5-10 GHz. It is interesting to note that the flux density below v a has 
the form F v oc v 2 , not the 5/2 spectral slope usually seen toward most radio 
sources. This is because the relativistic shock accelerates electrons to a power- 
law distribution (with energy index p given by N(7 e ) oc 7,T P ) above a minimum 
energy 7 m , which initially radiate their energy most of their energy at v m » 
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v a . The flux below v a depends only the angular size of the source and the 
fraction of the shock energy that goes into accelerating electrons pQ, and thus 
it is a useful diagnostic of the ratio of the energy of the shock and the density 
of the circumburst medium [E/n). 

As this optically thick radio source expands, a monotonic rise in the flux 
would be expected. It was therefore a considerable surprise when early obser- 
vations of GRB 970508 showed erratic, short term (~ hrs) and narrow band 
(~GHz) fluctuations in the centimeter emission ^j- The origin of these varia- 
tions |2] was traced to the scattering of the radio emission, owing to the small 
angular size of the fireball, as it propagates through the turbulent ionized 
gas of our Galaxy. This is a large and complex subject ^3E|> but for the 
purpose of this review it is sufficient to note that for typical lines of sight the 
modulation of the flux densities is near a maximum at frequencies near 5-10 
GHz. Coincidently, this is the same frequency range where v a typically lies 
and where the majority of radio observations are being made. While interstel- 
lar scintillation adds a certain degree of complexity to interpreting afterglow 
light curves, it also allows us to use the Galaxy as a large lens to effectively 
resolve the fireball. The observed "quenching" of diffractive scintillation from 
GRB 970708 four weeks after the burst |18II19| lead to estimate of the angular 
size, demonstrating superluminal expansion and providing an early confirma- 
tion of the fireball model. 

In many instances (201 EH E21 EI3 El bright, short-lived radio "flares" 
are detected at early times (t < 3d). The emission is much brighter than 
expected from a backward extrapolation of the light curve, and the level of 
fluctuation is too great to be accounted for by interstellar scintillation. One 
of the best-known examples is the radio flare of GRB 990123 [201 ^ which was 
accompanied by a 9th magnitude optical flash [251 - This prompt optical and 
radio emission is thought to be produced in a strong reverse shock which 
adiabatically cools as it expands back through the relativistic ejecta (2SJ- The 
strength and lifetime of this reverse shock emission is sensitive to the initial 
Lorentz factor, r o , of the shock and the density structure of the circumburst 
medium |271 |2"I] . To properly constrain these values requires that the peak 
of the emission be measured. This is difficult to do with optical observations, 
which require a response time on the order of the burst duration, while radio 
observations require a response time of only 12-48 hrs. 

On a timescale of days to weeks after the burst, the subsequent evolution 
of the radio afterglow (Fig. EJ can be described by a slow rise to maximum, 
followed by a power-law decay. The radio peak is often accompanied by a 
sharp break in the optical (or X-ray) light curves 28 , 22| ■ The most commonly 
accepted (but not universal) explanation for these achromatic breaks is that 
GRB outflows are collimated. The change in spectral slope, a, where F v oc 
t a yfi occurs when the r of the shock drops below 0J 1 , the inverse opening 
angle of the jet (301 E]. Since the radio emission at vr initially lies below the 
synchrotron peak frequency v m the jet break signature is distinctly different 
than that at optical and X-ray wavelengths. Prior to the passage of v m the 
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jet break is expected to give rise to a shallow decay t -1 / 3 or plateau t°, in 
the optical thin (y a < vr) or thick (y a > regimes, respectively. Another 
recognizable radio signature of a jet-like geometry is the "peak flux cascade" , 
in which successively smaller frequencies reach lower peak fluxes (i.e. F m oc 

1/2 

Vm )• Taken together, these observational signatures can be used to infer the 
opening angles dj of wide angle jets. Such jets are hard to detect at optical 
wavelengths because the break is masked by the host galaxy, which typically 
dominates the light curve between a week and a month after the burst |32II33| . 
Once the real geometry of the outflow is known [341 13*5] the energy released 
in the GRB phase and the afterglow phase can be determined. 

As noted above, the radio band is fortuitously located close to v a and 
as such it is a sensitive probe of the density structure of the circumburst 
medium. Extensive broadband modeling 36 has yielded densities in the range 
0.1 cm" 3 < n< 100 cm , with a canonical value of order n~10 cm . Such 
densities are found in the diffuse interstellar clouds of our Galaxy, commonly 
associated with star-forming regions. A density of order 5-30 cm -3 is also 
characteristic of the interclump medium of molecular clouds, as inferred from 
observations of supernova remnants in our Galaxy (e.g., Chevalier 1999 and 
references therein). Based on X-ray and optical observations alone, there have 
been claims of high n> 10 4 cm -3 """"3 EH] or low n<g; 10~ 3 cm~ 3 gH| cir- 
cumburst densities. However, in several of these cases when the radio data 
has been added to the broadband modeling (i.e. constraining v a ), there is no 
longer any support for either extreme of density |22l l3*3*j . 

One unsolved problem on the structure of the circumburst environment 
is the absence of an unambiguous signature of mass loss from the presumed 
massive progenitor star in afterglow light curves 0T] . Although there are some 
notable exceptions (e.g., Price et al. 2002), most GRB light curves are best 
fit by a jet expanding into a constant density medium instead of a radial 
density gradient, p oc r -2 |3fi| . Part of the solution may lie in reduced mass 
loss rates due to metalicity effects, or the motion of the star through a dense 
molecular cloud 03| , both of which act to shrink the radius that the pre-burst 
wind is freely expanding. It is equally likely that our failure to distinguish 
between different models of the circumburst medium is due to the lack of 
early afterglow flux measurements, especially at millimeter and submillimeter 
wavelengths where the largest differences arise [441 14"5] . The resolution of this 
conflict is important as it goes to the heart of the GRB progenitor question. 

At sufficiently late times, when the rest mass energy swept up by the 
expanding shock becomes comparable to the initial kinetic energy of the ejecta 
(~100 days), the expanding shock may slow to non-relativistic speeds 02]. A 
change in the temporal slope is expected at this time (Fig. EJ with onr — 
(21 — 15p)/10 for a constant density medium, independent of geometry. This 
dynamical transition provides a simple and power method to derive the kinetic 
energy of the outflow which has expanded to be quasi-spherical at this time. 
In contrast, most energy estimates made at early times require knowledge 
of the geometry of the outflow (401 1471 l3*B] . Using the late-time radio light 
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curves and the robust Taylor-Sedov formulation for the dynamics we can infer 
quantities such as the kinetic energy, ambient density, magnetic field strength, 
and the size of the fireball. The radius can be checked for consistency with the 
equipartition radius and the interstellar scintillation radius. This method has 
been used for GRB 970508 [IHJ and for GRB 980703 (Berger, priv. comm.), 
yielding energies of order few x 10 50 erg, in agreement with other estimates. 

Finally, the radio light curves at late times may flatten due to the presence 
of an underlying host galaxy. Most GRBs studied to date have optical/NIR 
hosts but only about 20% have been seen at centimeter and submillimeter 
wavelengths 0311111301. This radio emission, if produced by star formation, 
implies star formation rates SRF~ 500 M Q yr _1 and Lm > 10 12 L©, clearly 
identifies these GRB hosts as ultraluminous starburst galaxies which are all 
but obscured by dust at optical wavelengths. This is an emerging area with 
great potential for studying cosmic star formation with a sample of galax- 
ies selected quite differently than other methods. Preliminary studies have 
already shown that GRB-selected galaxies are significantly bluer than other 
radio-selected samples |50) . 
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